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Abstract 

We prove the existence of an ordered low temperature phase in a 
model of soft-self-avoiding closed random surfaces on a cubic lattice by 
a suitable extension of Peierls contour method. The statistical weight of 
each surface configuration depends only on the mean extrinsic curvature 
and on an interaction term arising when two surfaces touch each other 
along some contour. The model was introduced by F.J. Wegner and G.K. 
Savvidy as a lattice version of the gonihedric string, which is an action 
for triangulated random surfaces. 
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1 Introduction 



The gonihedric string was introduced by Savvidy et. al. ^ as a new 

action for random surfaces. For triangulated surfaces, the action reads 

^ = ^ E \x^-x,\e{a,,), (1) 

<i3> 

where 

e{aij) = \tt - aij\ (2) 

and Uij is the angle between the neighbouring triangles with common link 
< ij >. This string model can be considered as a natural extension of the 
Feynman integral over paths to an integral over surfaces in the sense, that 
both amplitudes coincide in cases, when the surface degenerates into a single 
particle world line. The simulation of (|l|) shows flat surfaces [Q], although 
some problems arise from the failure to suppress the wanderings of vertices 

i- 

One possibility to regularize the gonihedric action is to formulate the model 
on the euclidean lattice, which has been done by Wegner and Savvidy. There 
are two essentially distinct cases which correspond to non-self-avoiding surfaces 
[0] and to soft-self- avoiding ones |p. In the former case self- inter" -sections of 
the surface do not produce any additional energy, while in case of soft-self- 
avoiding surfaces, the energy contribution of a link, where self-intersection 
occurs was defined by times the number of pairs of plaquettes which meet 
under a right angle. It was shown, that in both cases a Z2 Ising model can be 
constructed, which reproduce the same surface dynamics and an extension to 
{d — n)-dimensional hypersurfaces on a d-dimensional lattice was also given. 

In this paper we consider the case of {d — l)-dimensional soft-self-avoiding 
surfaces on a d-dimensional cubic lattice. The statistical weight of each surface 
configuration is given hy E = I2 + 4:1^, where I2 is the number of links, where 
two {d — l)-dimensional plaquettes meet under a right angle (i.e. the mean 
extrinsic curvature) and Z4 is the number of plaquettes, where four plaquettes 
meet. This additional term arises, when two surfaces touch each other and 
is responsible for the soft-self-avoidance. The equivalent spin systems con- 
tain just ferromagnetic nearest neighbour and antiferromagnetic next nearest 
neighbour couplings. If one wants to allow arbitrary self-intersection coupling 
A; > 0, i.e. E = I2 + Akl/^, the corresponding spin-hamiltonian contains also 
a plaquette term. The two dimensional model for k = 1 does not seem to 
show a phase phase transition |8|, whereas in three dimensions a second 
order phase transition occurs at Pc ~ 0.44 which is close two the critical tem- 
perature of the two dimensional Ising model. These results were obtained by 
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numerical simulations. The three dimensional gonihedric Ising system was also 
discussed for arbitrary A; > in Q using mean field methods and simulations. 
The authors find a second order phase transition for k > and a qualitatively 
different behaviour for the case k = 0. The critical temperature for A; > was 
shown to increase with increasing k. 

In this work we prove the existence of an ordered low temperature phase 
for k > in d > 3 dimensions using a suitable extension of Peierls contour 



method |11]. After recapitulating the model, we adapt Peierls argument to 
the case k = 1. An upper bound for /3c for the three dimensional model is 
found, which is in agreement with [^, ^. Finally we show, how to generalize 
the argument for arbitrary k > 0. 

2 The Model 

Consider a d-dimensional euclidean lattice, with lattice points r G Z'^. We 
define a {d — n)-dimensional hyperplaquette J^ai,...,a„(r) (all Oj different) by 

^ai,...,a„ (r) 

= |x G i?'^ I = Tai, Ta < Xa < Ta + 1 for all a 7^ ajj (3) 

For convenience we call a (d — l)-dimensional hyperplaquette simply a pla- 
quette, a {d — 2)-dimensional hyperplaquette a link and a {d — 3)-dimensional 
hyperplaquette a vertex. 

A closed surface on the lattice is a collection of plaquettes, where at each 
link an even number of plaquettes meet. In a hamiltonian for closed surfaces 
M on the lattice is defined in the following way: 

Attach plaquette variables Up to each plaquette P of the lattice and define 

[ -1 if PGM 
^^ = 1 +1 if P^M 

The energy of M is given as a sum over all links 

H= J2 ^link, (5) 
all links 

where -ffiink contributes 4J if four plaquettes meet at the corresponding link, 
J if two plaquettes meet perpendicular and zero in all other cases. Note that 
4 is always the maximum number of plaquettes, which can meet at a link. 
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independent of the dimension of the lattice. In terms of the hnk variables this 
can be written as 



^link = \j{2 -Ui- U.i){2 -U2- U.2). 

U2 



I ^ - - 



U. 



(6) 



C/-2 



Ui, U-i, U2, U-2 can not be chosen independently since the surface is closed. 
This condition requires 

U1U-1U2U-2 = 1. (7) 

To resolve this constraint, spin variables a = ±1 were attached to the sites of 
the dual lattice. 

U2 




The plaquette variables can then be represented as 

?7l= (71(72, ?7_l = (7_i(7_2, U2=(yi(T-2, i7_2 = (7-1(72, 

and the full hamiltonian becomes equivalent to 
'd, 



H = jy: 



a<l3 



(8) 



(9) 
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In this expression r runs over all sites of the dual lattice, where the spins 
are located and 6^,6^ denote unit vectors parallel to the d cubic axes. One 
sees that the equivalent spin system contains ferromagnetic nearest neigh- 
bour and antiferromagnetic next nearest neighbour couplings, where the ratio 
between them is fixed and given by 

Jferromagnct _2(d-l). (10) 



J antiferromagnet 



For arbitrary self- intersection coupling k > 0, the energy contribution of a 
link, where four plaquettes meet is equal to AJk. In this case equation 
must be replaced by 

i^link = ^Jk{2-Ui-U^i){2-U2-U^2) 

+^Jil - k){2 - UiU-1 - U2U-2) (11) 
and the corresponding spin-hamiltonian becomes 



^{d - 1)^ - (d - l)fc ^ a{r)a{v - e,) 



a 



(a{r - e„)(7(r - e^) + <T(r)a(r - eo, - 



^ a<l3 



1 — k / 

— ^ [^(^)^(^ ~ ea)o-(r - e„ - e;3)cr(r - e/3 

a<f3 



(12) 



This hamiltonian also contains a plaquette term as long as A; 7^ 1. 

Because in (^) the plaquette variables were represented as products of 
nearest neighbour spins, we can restore the original surface M from the spin 
configuration by simply choosing all plaquettes between spins of opposite sign. 
In this sense, the surface M can be considered as a domain wall, which sepa- 
rates spins of positive from spins of negative sign. In contrast to the ordinary 
ferromagnetic Ising model, the energy of a spin configuration is not given by 
the total number of plaquettes of M, i.e. the surface, but hy E = {I2 + 4A: U)J, 
where I2 is the number of links contained in M, where two plaquettes meet 
perpendicular and I4 is the number of links, where 4 plaquettes meet. Be- 
cause of this, one can always insert plain domain walls by swaping a whole 
{{d — l)-dimensional) spin layer without changing the total energy, provided 
the inserted walls do not cross any existing surface. In particular the ground 
state of the finite lattice system containing A^'^ spins (d-dimensional cubus) is 
{d (2^ - 2) + 2)-fold degenerate. 
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3 Existence of a phase transition for d > 3 

In this section we will prove, that the model defined in the previous section 
shows a phase transition, provided d > 3. This will be done with a method 
which is due to Peierls The argument proceeds in two steps, (i) We 

show that the spontaneous magnetization M^v for a finite system with special 
boundary conditions is bounded from below, if the temperature is lower than 
a critical temperature T^, i.e. 

Mtv > a > for T < T^. (13) 

The hat indicates the special boundary conditions and the index refers to 
the finite system containing spins, a is independent of A^, but sensitive to 
the special boundary conditions we impose, (ii) The free energy /jv(T, //) is 
concave as a function of the magnetic field H, i.e. 

- fN{T,H)-UT,^) ^ ^^^^^ = 0) > « > 0, (14) 
H 

for T < Tc and H > 0. The thermodynamic limit liuiN^oo InIT, H) exists 
in the sense of van Hove and is independent of the boundary conditions. The 
limiting free energy is concave as well. This implies that equation (|T^ ) remains 
valid in the thermodynamic limit. Since the magnetization is antisymmetric 
in H, it follows that f{T, H) is not analytic m. H = Q \i T < T^- 

To get an idea how equation ( p!3[) can be shown, consider a 3-dimensional 
finite volume with spins and fix the spins at the boundary to +1. The 
ground state of this system is given by the spin configuration, in which all 
spins are positive, since then no domain walls are present and hence the total 
energy is zero. If we now swap a little 'island' of spins inside the volume to 
-1, we need an amount of energy which is essential proportional to the 
total number of edge elements of the domain wall, which was established by 
swaping the spins. Therefore, the probability of occurrence of a domain wall 
with I edge elements is of order v{l)e~^'^\ where v{l) is the total number of 
possible edge configurations with I edge elements. Now for d = 3, the total 
number of negative spins inside a domain wall is bounded by ( j^)'^ (volume of 
a cube), hence the mean number of negative spins will be bounded by 

(^-)^^E(^)'KOe-^^^ (15) 

where C is some constant. If does not increase too rapidly for increasing 
I, the right hand side of (|l5| ) will be smaller than ^ for large enough (3 > (5c- 
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Thus the spontaneous magnetisation will be strictly non zero, if /3 > (3c- The 

argument can easily be extended to d > 3, if (j^) replaced by { 2d(d-i) ) ''"^ • 
It does not work in two dimensions, since in this case, the edge elements are 
vertices and the total number of negative spins inside a domain wall with a 
given number of vertices can be arbitrary large. Indeed, in two dimensions, 
the model does not seem to show a phase transition We now proceed 

with the details for the case k = 1. 

Consider a finite volume of the dual lattice containing N spins in its interior 
V \ dV and fix the spins at the boundary: 

cj(r) = +1 if redV (16) 

If a certain configuration cr = (fji, . . . , cjtv) of the spins is given we denote by 
M[cr] the corresponding closed surface on the original lattice which separates 
regions of negative from regions of positive spins. We call a link contained in 
iVf[<T] an edge element, if one of the following two conditions is fulfilled: 

• Four plaquettes of M[cr] meet at the given link. 

• Two plaquettes of M[cr] meet perpendicular at the given link. 

The set C*°* of all edge elements of M[<t] can be divided into connected edge 
diagrams Ci , . . . , in a unique way 

C*°* = CiUC2U ...UCn (17) 

The energy E[cr] of the spin configuration can then easily be expressed in terms 
of connected edge diagrams since M[cr] is the original interacting surface and 
therefore only the edge elements of M[cr] contribute to E[cr]. We define the 
energy contribution [C] of a connected edge diagram C as 

E[C] = Jih + ^h), (18) 

where I2 denotes the number of edge elements of C, where two plaquettes meet 
and U the number of edge elements, where four plaquettes meet. The total 
energy of a spin configuration can then be written as 

E[a]=J2E[C]. (19) 

c 

Clearly, the number of all connected edge diagrams with a given energy con- 
tribution E that fit into the finite volume V is finite. We denote this number 



6 



by g{E). Next we attach variables xfl""];* = ■ ■ ■ idiE) to all these con- 
nected edge diagrams Cf. xfW] assume the value +1 if the corresponding 
edge diagram appears in a and otherwise. 

The partition function for the spins is given by 

Zn = Y.^-^'^^'"^ (20) 
a 

We are interested in the thermodynamical expectation value {xf^ that is 
the probability of appearance of Cf: 

(xf),=^i:xf[.i.-^^"" (21) 

In this expression only spin configurations which contain Cf" contribute. Num- 
bering them by <ti, . . . , a-^^ (|2l|) reads: 



Consider a particular spin configuration aj in which Cf appears. The con- 
nected edge diagram Cf belongs to a domain wall, which encloses a certain 
number of spins. By reversing all these spins, the domain wall disappears and 
the resulting configuration <t| will not contain Cf. Therefore 

E[aj]=E[a*] + E[Ch (23) 
From this we obtain the following estimate: 

Zn>Y. e^'^^'t'^^*' > e^^I^^^^l ^ e-^^['^^l . (24) 
j=i j=i 

Together with (^) this implies: 

(xf)^<e-^^[^^"l. (25) 

We can now find an upper bound for the mean number of negative spins inside 
V as follows: 

Consider a given configuration cr and denote its connected edge diagrams 
by Ci , . . . , Cn • Each diagram belongs to a domain wall, which contains a 
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number rij of negative spins. If Ij < E[Cj]/J denotes the number of edge 
elements in Cj, rij is bounded by 

Ij \^ f E[Cj] \^ 



This relation becomes an equation, if the domain wall is a simple cube. The 
total number of negative spins [a] is therefore bounded by 

Using the variables xfi""]) we can write 

9(E) , „ s 

~-WSi:ExfW(55(^) , (28) 

where the first summation is over all possible energies E, a connected edge dia- 
gram can assume. From ( ^5|) we get an upper bound for the thermodynamical 
expectation value of N-[cr]: 



Note that the right hand side of (|29|) is not defined for d = 2. To proceed 
with the argument, we need an upper bound for g{E). Consider the following 
construction method, which can be used to build up every possible diagram: 

1. First we number all the vertices {{d — 3)-dimensional hyperplaquettes) 
of the lattice and keep this numbering once for all fixed. 

2. Next we choose one of the N lattice points and attach 3 edge elements, 
pointing in three given directions. For example if p = (xi, . . . ,Xd) de- 
notes the chosen point, we can attach the edge elements 

Ci = (xi,X2,A3,...,Arf), 

C2 = (xi,A2,X3,A4,...,Arf), 

C3 = (Ai,X2,X3, A4, . . . , Arf), 

Xi < Aj < Xi + 1. This can be done since such an edge configuration 
must occur at least once in any possible connected edge diagram. We 
think of this vertex as a corner of the closed surface. This fixes the 8 
spins surrounding the vertex p. 
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3. Imagine that a connected subdiagram exists already. Now we choose the 
vertex with the lowest number, which occurs in the construction so far, 
i.e. where at least one edge element ends and where the configuration 
of the edge elements surrounding this vertex has not been specified yet. 
Next we do this specification by choosing the values of not more than 
four of the eight spins, surrounding the vertex, since we know that at 
least one edge element ends at this vertex already. 

4. Finally by repeating construction step 3., all possible connected edge 
diagrams can be constructed. 

To find an upper bound for g{E), we estimate the maximum number of pos- 
sible outcomes of the construction procedure just given. If the considered 
finite lattice volume contains N spins, we have N choices for construction 
step 2. Each time we perform construction step 3, we have different numbers 
of choices how to specify the remaining spins. This depends on how many 
edge elements and what kind of edge elements are already ending at the ver- 
tex under consideration. The maximum number of choices arises, if just one 
edge element, surrounded by two or four plaquettes, is already present at the 
vertex, leaving four spins to specify, i.e. 16 possibilities. These possibilities 
are shown in figure 1, together with the corresponding energy contribution of 
the specified vertex. Edge elements where two plaquettes meet are indicated 
by a simple edge. Edge elements, where four plaquettes meet are indicated as 
double edges. The left column shows all cases, where one simple edge coming 
from the left was already present and correspondingly the right column shows 
all cases, where one double edge coming from the left was already present. 
Prom this picture we read of the maximum numbers of choices n{E), how to 
specify the spins, such that a vertex of energy E arises: 



maximum number 
of choice n{E) 




1 




4 


X^ 


2 


X^ 


2 


x^ 


4 


x' 


3 


x^ 


4 


x^o 


4 




1 
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The function Gn{x) = Nx^ f^{x), where /(x) is defined by 
fix) = Y.<E)x'"' 

E 

= x2 + 4x3 + 2x'' + 2x^ + 4x^ + 3x^ + 4x^° + 4x^^ + x24 (30) 

can now be interpreted as a generating function which counts the number 
of all connected edge diagrams with energy E at least once, which can be 
constructed by repeating construction step 3 n-times, i.e. the coefficient of 
^2{d-2)E/j expansion of Gn{x) is an upper bound for the number of 

those diagrams. Since ah edge elements were counted 2{d — 2) times, we have 
to calculate the coefficient of x^^'^~'^^ ^/'^ and not simply x^/'^ . Summing Gn{x) 
over all possible numbers of steps 

00 3 
G{x) = Y.Gn{x)=N-^—- (31) 

n=0 ^ J^^> 

we get a generation function, which counts all connected edge diagrams at 
least once. Therefore 



^ Q2(d-2)E/J 

^^^^ - {2{d-2)E/J)\ 9x2{'^-2)i?/j '^(^^ 
provides an upper bound for g{E). We can write G{x) as 



(32) 



x=0 



24 

G{x) = Ny (33) 



(34) 



Putting this expression in equation (|3^), we obtain 

\E/J 

-2) 



k=l ^ \Xu 



If we denote the smallest \xy\ by Xmim which clearly dominates the increase 
of g{E) for large E, we can further approximate this expression as follows: 



2) 

mm 



E/J 

I =:iVac('^-2)i^/^ (35) 



where 



a < 0.625 
c < 3.882 
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Substituting this result into equation (29), we find an upper bound for the 
density of negative spins: 



N -^\2d{d-2)j) '''^ ^^^> 

This inequahty remains vahd in the thermodynamical hmit. The sum on the 
right hand side converges, if /3 is larger than ^j^ln(c) and approaches zero, 
if (3 goes to infinity. Therefore for /3 > /3c, it will be smaller than ^, which 
implies a non zero spontaneous magnetisation. This proves the existence of a 
phase transition. 

Equation (p9|) together with an upper bound for g{E) given in (|^) also 
provides a lower bound for the critical temperature Tc- We find 

Tc > 0.67 J for d = 3 (37) 
Tc > 0.36 J for d = 4 (38) 

The value for (i = 3 is in agreement with the critical temperature found in 
|, |], which was Tc « 2.28 J 

The method used above to prove the existence of an ordered low temper- 
ature phase for k = 1 can easily be generalized to arbitrary A; = | > 0, where 
p,q & N. We have 

E 1 

— = l2 + AkU>l = l2 + li if k>- (39) 

ifj = | + i4>< = fa + <. if k<\ (40) 



Therefore equation (29) remains unchanged if A; > whereas in case k < ^ 
we have to replace it by 

An upper bound for g{E) can be found by expanding the generating function 

G(x) = iV^, (42) 



+ 4x^ + 2x* + 2x^ + 4x^+''9 + 2x 



4+4£ 



+x 9 + 4x 9 + 4x 



122i , 242 
9 + X 9 



(43) 
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in powers oi xi , i.e. 

GM=^E^=£fwi:(-s)(i)"].? (44) 

j XI -Xj n=0\ j -^j' ^-^j' J 

Prom this we obtain 

The estimate for g{E) becomes smaller for increasing k and diverges for k ^ 0. 
For A; ^ oo the behaviour of g{E) for large E is dominated by the smallest 
zero of the function 1 — x^ — Ax^ — 2x^ — 2x^, which gives 

g{E) < const {3.74)^"^-^^^/-^ for A; ^ oo (46) 



4 Summary 

We have proven the existence of an ordered low temperature phase for the 
spin systems defined in section 2, which can be considered as the lattice reg- 
ularization of the gonihedric string. Each spin configuration corresponds to a 
particular random surface which is simply the domain wall separating regions 
of spin with opposite sign. This observation was crucial, since it allowed us to 
apply Peierls contour method. However the naive application of the argument 
would yield an inequality of the following form in d = 3 dimensions: 

3 

(iV_)^ <12{iy ^(^) e-^'^rmnib)^ (47) 

where {N-)j^ is the expectation value of the number of negative spins, calcu- 
lated for the finite system with boundary condition cr(r) = +l,r G dV. The 

3 

factor (6/6)2 is the maximum number of negative spins that a closed surface 
with b plaqucttcs can enclose and i^{b) denotes the number of closed surfaces 
with b plaquettes. This number is known to grow exponentially, in fact 

u(b)<N3''-\ (48) 

For Emini^) we have to put the minimum energy of a closed domain wall with 
b plaquettes which is 

E^Ub) = l2\r- {d = 3) (49) 
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This shows, that the right hand side of (|4^) actually diverges for every (3 due 
to the exponential growth of i^{b). We therefore had to modify the argument 
as discussed in section 3. This was possible since the energy of a closed surface 
is essentially proportional to the number I of edges and the number of negative 
spins inside is bounded by (j^)^- Instead of ( p7| ) we arrived at 

for < < I and 

(51) 

for k > J. Finally we proved that the number of connected edge diagrams 
with given energy g{E) does not grow faster than exponentially for any k > 0. 

As pointed out before, the argument does not work in two dimensions, 
because in this case the (0-dimensional) edges are not connected and therefore 
do not give a restriction for the number of negative spins. 
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